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ABSTRACT 

Elasticity of the liner and its effects on cavitation were numerically 
determined for a squeeze film damper subjected to dynamic loading. The 
loading was manifested as a prescribed motion of the rotof undergoing 
noncentered circular whirl. The boundary conditions were implemented using 
Elrod's algorithm which conserves lineal mass flux through the moving 
cavitation bubble as well as the oil film region of the damper. Computational 
movies were used to analyze the rapidly changing pressures and vapor bubble 
dynamics throughout the dynamic cycle for various flexibilities in the damper 
liner. The effects of liner elasticity on cavitation were only noticeable for 
the intermediate and high values of viscosity used in this study. 

INTRODUCTION 

Earlier work by Jacobson and Hamrock (Refs. 1 and 2) using high-speed 
photography revealed the existence of vapor cavitation in dynamically loaded 
journal bearings. A polymethylmethacrylate (PMMA) housing was used for 
visualization of the cavitation. In their studies they investigated the 
effect of roller material, the influence of L/D, effect of static and dynamic 
eccentricities and influence of speed. A theoretical model implementing the 
Elrod algorithm (Ref. 3) was developed by Brewe (Ref. 4) for the dynamically 
loaded journal bearing. The algorithm conserves mass throughout the 
computational domain and essentially mimics the Jakobsson-Floberg and Olsson 
boundary conditions (Refs. 5 and 6). Later some of the earlier experiments 


were used to corroborate the theoretical model (Ref. 7). A prediction of 
32 ms bubble life for nonstationary cavitation was obtained from the 
theoretical model. Using high-speed photography for the same set of 
conditions, a bubble life of 23 ms was observed in the laboratory. This 
discrepancy was attributed to deformation and temperature effects as well as 
some uncertainty in correlating actual life with duration of visibility. To 
achieve a more representative comparison it is necessary to Include the 
thermal effects in the model as well as deformation of the shaft and housing. 
This analysis represents a first step toward including first-order deformation 
effects with the Elrod algorithm. As such the model assumes that the bearing 
is essentially a thin tube surrounded by a rigid housing. Since the periphery 
of the bearing is much larger than its thickness, the radial deformations of 
the latter at a point may be assumed to be proportional to the pressure at the 
point. This model is commonly referred to as the Winkler foundation and has 
been used by numerous investigators, i.e., (Refs. 8 to 12). At approximately 
the same time that Dowson and Taylor (Ref. 9) published their work for a 
compliant thrust bearing application, Castelli, et al . (Ref. 13) published 
their findings using an exact analytical approach. In this paper and in a 
discussion to the Dowson paper (Ref. 9), the limitations of the Winkler 
foundation are discussed. Based on this, the thin liner materials used in 
this investigation appear to be suited to the assumptions of the theory. 
Realistic edge effects are not expected however. 

The purpose of this paper is to study the elastic effects of the bearing 
liner on cavitation in a squeeze film damper. It has been shown for steadily 
loaded bearings that increasing the liner flexibility increases the cavitation 
angle (Refs. 14 to 17) using Reynolds boundary conditions. It is believed 
that under dynamic conditions the elasticity of the liner may influence the 
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bubble dynamics and hence the damping characteristics of the squeeze film 
damper. 

NOMENCLATURE 


B parameter, 


6nw^ d 


1 - 


2 a 


1 - 


(AR/R )ER V 
D shaft diameter, m 

d thickness of the bearing liner, m 
E Youngs Modulus, N/m 2 

e eccentricity, m 

g switch function, cavitation Index 

H dimensionless film thickness, h/AR 

h film thickness, m 

L length of the bearing, m 


L/D length to diameter ratio 


lineal mass flux, kg nr 1 s 


-1 c -1 


p fluid pressure, N m -2 

p a ambient pressure, N m~ 2 

p c cavitation pressure, N m -2 

. , (AR/R) 2 p 

p dimensionless pressure, ~ 

R radius of shaft, m 

t time, s 

U sum of the surface velocities In x direction, m s _1 

V sum of the surface velocity vectors, m s - ^ 

W load capacity, N 


w dimensionless load capacity. 


(AR/Rr M 
6nuj 
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liquid bulk modulus, N nr 2 15 
Y angular position of minimum film, rad 

A dimensionless radial deflection', S/AR 

AR radial clearance, m 

6 radial deflection, m 

e eccentricity ratio, e/AR 

0 fractional film content in cavitation zone; density ratio (p/p c ) in full 
film zone 

p dynamic viscosity, N s nr 2 

p fluid density within full film zone, kg nr 2 

p c fluid density within cavitated zone, kg m~3 

<p angular coordinate relative to minimum film line, rad 

uq orbital angular velocity of journal center about fixed point relative to 

housing center, rad s -1 

w s angular velocity of journal about its own center, rad s -1 

THEORY 

Cavitation Algorithm 

A computer program was developed based on Elrod's algorithm (Ref. 3) that 
combines a control volume approach to deriving the finite-difference equation 
with a switching function that automatically handles the moving boundary 
condition at the liquid-vapor interface. The complication in satisfying the 
boundary condition is avoided by defining a single independent parameter 9 
which links the full-film region and the cavitated region. This parameter is 
defined by the ratio of the local density to that of the cavitation density, 

1 .e. , 
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It follows that in the cavitation region 0 < 1.0 and In the full-film 
region 9 > 1.0. A switch function g was defined as "on" (g = 1) and "off" 

(g = 0) in the full-film and cavitation regions respectively. By using this 
property, a universal Reynolds equation was defined as Indicated: 

+ $ • iH = o (1) 

where m represents the lineal mass flux and is given by 

^ 6gce>?r (2) 

P c 2 12 M 

In the cavitation region, the second term vanishes since g = 0. This 
equation was Integrated in time by using an Euler's method to advance in 
time. At every time step the difference equation was solved by an Alternating 
Direction Implicit (ADI) scheme. The switch function was updated after each 
half-time step. Occasionally, during the collapse of the vapor bubble, liquid 
was piling up at the boundary. This was due to an Inherent time lag of the 
switch function with the calculation. Consequently, the switch function was 
immediately updated to accommodate the movement of the boundary, and the 
calculation was reiterated. 

The results of the program developed were authenticated by comparison to 
an experiment of Ref. 6. When compared to the pseudo-Giimbel boundary 
conditions (Ref. 4), the results predicted a significant difference in the 
computed load capacity (up to 20 percent). Furthermore, it was found that the 
size, life, and location of the cavitation bubble were considerably affected 
by imposing the pseudo-GLimbel boundary conditions. 

Fluid Velocity Under Dynamic Loading 

Consider the bearing geometry shown in Fig. 1. As shown, the minimum-film 
thickness rotates clockwise around the bearing at a whirl frequency of y. 
Jacobson and Hamrock (Ref. 2) show that from this diagram y can be evaluated 
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(3) 


tan 


1 e^j sln(u^t) 

e„ + e . cos(u.t) 
s a d 


so that 


e d Ec s co,< “d t> + e d 3 “d 
Y * 2 2 

c + e . + 2e„e . COS(w.t) 
s d s d d 


(4) 


The sliding velocity of the journal is 

U = (wj - 2 y)R 
Film Thickness and Deformation 

The film thickness for a bearing having a rigid liner can be given as 

h = AR - e cos <p 

where the total eccentricity e defined as the instantaneous distance 
between Ob and Oj is 

e d s1n ( “d^ 
e * sin y 

If one considers deformation of the bearing surface, the film thickness 
can be written as: 


(5) 


( 6 ) 


(7) 


( 8 ) 


(9) 


h = AR - e cos <p + 6 

or normalizing with respect to the clearance, AR, we have 

H = 1 - c COS cp + A 
where S represents the radial deflection and A = S/AR. 

This study considers only 1st order effects of the elastic deflection of 
a soft bearing l|ner on damper performance. These deflections are calculated 
using the assumptions adopted by Higginson (Ref. 8) for a soft liner material 
rigidly enclosed, i .e . , 

A = Bp (10) 

This simple expression can be derived by assuming the liner to be an 
elastic material that deforms linearly due to local pressures under conditions 
of plane strain. Accordingly, the parameter B can be shown to be 
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( 11 ) 


B = 6na) d t L _ ll ' 

(AR/R) 3 ER V 1 ■ 

- (AR/R) 2 p 

P 61110^ 

Since the film thickness is directly proportional to the local pressure, 
the solution of Eq. (1) requires nonlinear techniques. An iteration procedure 
was introduced that required the sum of the pressure field to deviate by no 
more than 0.05 percent due to iterations in film thickness. 

GEOMETRY AND OPERATING CONDITIONS 

The damper (bearing) simulated in this study is based on an experimental 
study of a journal-bearing squeeze film damper submerged In Dexron lubricant 
(Ref. 7). The range of operating conditions are listed in Table I. 

RESULTS AND DISCUSSION 
Effect of Flexibility on Cavitation 

The effect of flexibility on vapour cavitation was studied for a journal 
undergoing noncentered circular whirl. Three liner materials were used in 
this investigation to achieve the range of flexibilities. These materials 
were steel, polymethylmethacrylate (PMMA) and polytetraf luorethyl ene (PTFE). 
Table II lists these materials and their elastic constants in order of 
increasing flexibility for a range of different viscosities. The percent of 
time that a vapor bubble is predicted for one full orbit of the journal is 
listed under "Bubble life." To the right of this column are listed bubble 
inception and collapse. These values are given in milliseconds of time 
relative to the moment that the minimum eccentricity of the orbit occurs. 

Thus at startup, e = 0.1 and t = 0.0. The orbital period of the journal is 
67.7 ms during which time the eccentricity ratio varies between 0.1 and 0.9. 
Note that as the viscosity increases the bubble life increases to the extent 


7 


that cavitation exists throughout the entire cycle for both the steel liner 
and the PTFE liner. However, the effect of increased flexibility retards the 
sensitivity of bubble life to Increases in viscosity. This appears to be 
attributed to the mechanism leading to the collapse of the bubble. From the 
table we see that for the steel liner the bubble collapse is delayed as 
viscosity is Increased. But for the PTFE liner the bubble collapse appears 
insensitive to viscosity changes. 

From the standpoint of a viscous damper, it is logical that cavitation 
can only serve to decrease the damping characteristics. The above results 
lead one to infer that as the viscosity is increased better damping 
characteristics can be achieved using a PTFE liner over a steel liner for a 
limited range of viscosities. Further, Increasing the viscosity should give 
better damping characteristics, but the accompanying Increased bubble life 
could deter that desirable effect. 

Cavitation and Pressure Distribution 

Figures 2 and 3 illustrate the responding pressure variation to the 
motion of the shaft undergoing noncentered circular whirl for two contrasting 
situations. Figure 2 represents the events taking place for a rigid (steel) 
liner with parameter B = 0.000038. Figure 3 contrasts those same events for 
a soft (PTFE) liner having B = 0.004440. In both cases, the viscosity was a 
constant 0.099 Pa • s. In each figure, part (a) Illustrates the starting 
position (i.e., e = 0.1). Parts (a) to (f) depict several key events in the 
clockwise circular motion of the journal center within the bearing. 

Figure 2(a) represents the position of the journal within the housing and 
the associated pressure distribution at the instant that the eccentricity is 
at a minimum. The dark oblong region Indicates the shape and position of a 
vapor bubble downstream of the minimum-film line. The position of the bubble 
within the small clearance space of the journal bearing is also Indicated. At 
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this particular instant the tendency of the pressure flow along the side of 
the bubble to cause the bubble to collapse becomes greater since it Is 
proportional to h^. The journal center moves In a clockwise orbit from 
Figs. 2(a) to (b) (nearly one-tenth of a revolution). While the eccentricity 
Is increasing, the minimum film Is decreasing which gives rise to the squeeze 
action in the vicinity of the minimum-film line. A slight rise In the 
pressure field can be observed there. Note In Fig. 2(b) the collapse of the 
bubble occurs ~6 ms later than in 3(b). 

Cavitation Is absent from Fig. 2(b) to (c). In Fig. 2(c), a new 
cavitation bubble appears, and the pressure generation in the region of the 
minimum-film line becomes significantly large because of the squeeze effect 
taking place there. Figure 2(d) represents the configuration of the bearing 
when the maximum pressure occurs. Note that this happens before the maximum 
eccentricity Is reached. 

Comparing Fig. 2(d) with Fig. 3(d) we see that the peak pressure occurs 
at essentially the same eccentricity (I.e., 0.88) but that the shape and 
magnitude of the pressures are markedly different. For the PTFE liner In 
Fig. 3(d) the peak pressure Is 51 percent lower than the rigid liner of 
Fig. 2(d). Also it is illustrating the effect of the displacement In the high 
pressure region producing a flattening shape along the axial direction. 

Figure 2(e) represents the end of the squeeze cycle and the beginning of 
the separation cycle. Now the minimum film begins to Increase which creates a 
suction effect, an accompanying growth in the bubble, and a dissipation of the 
pressure hump noted in Fig. 2(f). It is interesting to note that the 
eccentricity is the same (0.88) for Figs. 2(d) and (f). The drastic 
difference in the pressure plot is a result of the squeeze velocity having 
opposite signs in the two situations. The minimum-film thickness continues to 
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increase as it moves in a clockwise direction. The representation shown in 
Fig. 2(a) Is repeated as are the representations for each succeeding orbit. 

Dimensionless Load Capacity Versus Eccentricity 
Figure 4 represents the dimensionless load capacity (w) versus 
eccentricity ratio (e> for different B parameters. By presenting the 
results using dimensionless load, (i.e., duty parameter) the effect of 
viscosity is normalized out and one can focus on the relative elastic effects. 
The predominant effect to note in this figure is that as B increases the peak 
dimensionless load decreases. These curves were generated for two different 
liner materials (i.e., steel and PTFE). The steel liner (B = 0.000025) is 
represented with the curve having the maximum peak load. The other two curves 
were generated with PTFE as the liner material. The reduced dimensionless 
peak load for the PTFE material is a result of larger deflections due to the 
local pressures. Thus higher film thickness in the region of higher pressures 
result in a pressure relaxation and consequent reduction in load capacity. 

This is a commonly observed result. 

CONCLUSIONS 

A study of the effect of bearing liner elasticity on vapor cavitation and 
damper performance under dynamic conditions is presented. The motion 
considered in this study was that of a journal undergoing noncentered circular 
whirl. The bearing liners considered were steel, polymethylmethacrylate, and 
polytetraf luorethylene. The analysis was performed using the Elrod algorithm 
for an Isoviscous lubricant. A range of viscosities were investigated (i.e., 
0.066 Pa • s to 0.66 Pa • s) to help elucidate the effect of liner 
elasticity. The study revealed the following: 

1. For the low values of viscosity in this study there was no effect on 
cavitation due to increasing the flexibility by two orders of magnitude. As 
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viscosity was Increased the effect of flexibility became apparent by 
Inhibiting the Increase In bubble life. 

2. Increasing the viscosity led to an Increase In bubble life for all 
liner materials considered. 

3. As the viscosity Is Increased, better damping characteristics can be 
achieved using a PTFE liner over a steel liner for a limited range of 

vl scosl ties . 

4. Increasing the flexibility of the bearing liner results In a reduction 
and flattening of pressures along the axial length of the bearing. 

5. As the parameter B Increases the peak dimensionless load decreases. 
This parameter essentially normalizes out the effect of viscosity. 
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TABLE I. - OPERATING CONDITIONS 


AR , m 

S.OxlCT 4 

R, m .......... 

0.0425 

L/D ..... 


£ . 

0.1 to 0.9 

co^, rad/s ....... 

-92.7 

« s , rad/s 

0 

(J, N/m 2 

1.72xl0 9 

P a , N/m 2 

1. 0133x10 s 

P c . N/m 2 

....... 0.0 


TABLE II. - EFFECT OF FLEXIBILITY ON CAVITATION 


Li ner 
material 

Materi al 

constants 

Viscosity, 

n(PaS) 


Bubble 1 i fe, 
percent cycle 

Inception , 
ms 

Col lapse, 
ms 

a 

E(GPa) 

Steel 

0 

30 

PO 

o 

a> 

0 


0.025 

77.13 

22.74 

7.23 








.029 

81 .71 

21.74 

9.34 







.099 

.038 

90.55 

19.44 

13.04 







.264 

.100 

100.00 









.660 

.250 

100.00 

— 


PMHA 

0.40 

2 

76 



76.73 

22.80 

7.03 

PTFE 

0.46 

0 

50 

1 

ESS 

76.25 

23.04 

6.94 








3.454 

77.87 

21.74 

6.94 







.099 

4.440 

81.26 

19.64 

6.94 







.264 

11.841 

100.00 

— 

— 







,660 


100.00 

— 
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